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Using geometric methods for linearizing systems of second order cubically semi-linear 
ordinary differential equations, we extend to the third order by differentiating the 
second order equation. This yields criteria for linearizability of a class of third order 
^ ' semi-linear ordinary differential equations, which is distinct from the classes available 

\l , in the literature. Some examples are given and discussed. 

1 Introduction 

' The study of non-linear differential equations initially centred about procedures to reduce 

■ them to linear form. The earliest attempts tried to approximate the non-linear equation 

by a linear one and use the solution to iteratively improve the approximation [15]. One 
\ then has to hope that the series converges. It is never very clear in this procedure in 

^ ■ what regime the essence of non-linearity will be lost in approximation. Lie used his point 

^ ■ transformations [8] to determine the class of scalar second order semi-linear ordinary dif- 

ferential equations (odes) that would transform under it to linear form [9], thus providing 
exact solutions of the odes. Such odes are called linearizable in the sense of Lie. More 
generally, one allows other transformation of variables. We shall simply call them lineariz- 
able. He obtained specific criteria for linearizability and obtained the most general form of 
the equation to be cubically semi-linear, apart from obtaining an algebraic classification of 
the equations that could be linearized. Some work on classification of systems of two such 
equations was done [17] and it was found that there are five classes of such linearizable 
systems. Some further work proved that for the scalar third order there are three classes 
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[10]. However, till recently nothing had been obtained about explicit linearizability criteria 
of systems of second order, or of third order odes. 

Chern [2,3] used contact transformations to discuss the linearizability of scalar third 

order odes reducible to the linear forms u'" = and ii!" + u = 0. Point transformations 
were used by Grebot [5,6] but were restricted to the class of transformations t = (f){x), 
u = il>(t, x). This work was generalized by Neut and Petitot, and Ibragimov and Meleshko 
[16,7] to the third class known to exist [17], u'" + a{t)u = 0. They obtained a larger class 
of third order semi-linear odes that is linearizable. Meleshko [14] also considered arbitrary 
point transformations to reduce third order odes of the form y'" = F\y, y', y"] and used 
the d/dt symmetry to reduce the order to two. He then used the Lie linearizability criteria 
to determine the linearizability of the scalar third order equation. 

Using the connection between the isometry algebra and the symmetries of the system 
of geodesic equations [4], linearizability criteria were stated for a system of second order 
quadratically semi-linear odes, of a class that could be regarded as a system of geodesic 
equations [11]. The criteria came from requiring that the coefficients in the equations, 
regarded as Christoffel symbols, yield a zero curvature tensor. We use the projection pro- 
cedure of Aminova and Aminov [1], which appeals to the fact that the geodesic equations 
do not depend on the geodetic parameter, to reduce the dimension by one and convert 
from a quadratically semi-linear system to a more general cubically non-linear system [12]. 
When applied to a system of two dimensions we get a single cubically semi-linear ode. The 
linearizability criteria so obtained are exactly the Lie criteria! Applied to a system of three 
dimensions, one obtains a system of two cubically semi-linear odes with extended Lie crite- 
ria. The system of odes so obtained is in the class of maximally symmetric equations (out 
of the total of five classes mentioned earlier). Some other classes come from projection of 
the original space to lower dimensional spaces. 

In this paper we consider the linearizability criteria of a class of scalar third order 
odes obtained by differentiating the cubically semi-linear system of second order odes and 
stating criteria of linearizability for the class so obtained. (It is worth mentioning that we 
obtain the same criteria if we first differentiate the system and then project it.) The class 
obtained here is distinct from all the earlier classes considered as there is no guarantee 
that there will be three arbitrary constants appearing in the solution. 

The plan of the paper is as follows. In the next section we present the geometrical 
notation used and briefly review the essential results used in the sequel. In section 3 we 
derive the linearizability criteria for the scalar third order ode in the form involving the 
second derivative and in a form quintic in the first derivative. In the following section 
we present some illustrative examples. Finally, in section 5 we give a brief summary and 
discussion of the results. 



2 Notation and review 

The following are well-known and can be found in text books. We use the Einstein 
summation convention that repeated indices are summed over the entire range of the 
index. Thus, A"'Ba stands for S"^^y4"i?a. The metric tensor will be represented by the 
symmetric (non-singular) matrix gij and its inverse by g^^ . The Christoffel symbols are 
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given by 

— 2^*'"^'"''= ~'~ 9km,j — 9jk,m), (1) 

where ,k stands for partial derivative relative to x^, etc. The Christoffel symbols are 
symmetric in the lower pair of indices, F*^ = T\j. 



In this notation, the system of n geodesic equations is 

- jk 



x' + rkx^i^ = 0, i,j,k = l,...,n, (2) 



where is the derivative relative to the arc length parameter s defined by ds^ = gijdx^dxK 
The Riemann tensor is defined by 

^jkl — ^ jl,k ~ ^ jk,l + mk^ jl ~ ^ ml^ jk^ l-^j 

and has the properties 

^]kl = -^]lk^ (4) 

^]kl + ^klj + ^Ijk = 0; (5) 

and 

^]kl;m + ^jlm;k + ^]mk;l = 0- (6) 

The Riemann tensor in fully covariant form is 

Rijkl = diml^kli C^) 

and satisfies 

Rijkl = —Rjiki- (8) 

The linearization criteria of [11] are Rijki = 0. 

For a general system of two odes for two functions, if it can be regarded as a system of 
geodesic equations, 

x" = a{x, y)x'^ + 2b{x, y)x'y' + c{x, y)y'^, 

y" = d{x, y)x'^ + 2e{x, y)x'y' + f{x, y)y'^ (9) 
the Christoffel symbols in terms of these coefficients are 

rl^ = -d,Tl^ = -e,Tl^ = -f. (10) 
The linearizability criteria are 

tty — bx + be — cd = 0, by — Cx + {ac — b^) + {bf — ce) = 0, 

dy-Sx- {ae - bd) - {df - e^) = 0, (6 + f)x = (a + e)y, (11) 
with constraints on the metric coefficients: gu = p,gi2 = 521 = Qt922 = 

Px = -2{ap + dq), Qx = -bp - (a + e)q - dr, Vx = -2{bq + er), 

Py = -2{bp + eq),qy = -cp - {b + f)q - er, Vy = -2{cq + fr). (12) 
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The compatibihty of this set of six equations gives the above four hnearization conditions 
(11). For a system of three equations we get eighteen such equations out. In general 
there are n^(n + l)/2, n > 2. One now obtains the required hnearizing transformation 
by regarding the variables in which the equation is linearized as Cartesian, thus having 
9ii = 922 = 1 and gi2 = 921 = 0, and looking for the coordinate transformation yielding 
the original metric coefficients [11]. 

Following Aminova and Aminov [1] and projecting the system down by one dimension, 
the geodesic equations become [12] 

x" + AbcX^x'^x'' + B^^x^x" + C^x^ + = 0, a = 2, . . . , n, (13) 

where the dot now denotes differentiation with respect to the parameter x^ and the coef- 
ficients in terms of the F^^s are 

Ac = -r^c) -Sfec = F^j, - 2(5"^F^)^, 

= 2F?, - <5,«Fh, D« = F?i, a, 6, c = 2, . . . , n, (14) 

which is cubically semi-linear. For n = 2, writing h = a — 2e,g = f — 2b, we get the Lie 
criteria for the scalar equation involving the two auxiliary variables b and e, the former of 
which corresponds to our b and the latter to our — e: 

3bx — 2gx + hy — 3be — 3cd = 0,by — Cx + b^ + bg + ch — ce = 0, 

dy + ex — bd — dg — + eh = 0, 3ey — 2hy — gx + 36e -|- 3cd = 0, (15a) 
which yield, in terms of the four coefficients of the cubically semi-linear equation [12], 

3{ch)x + Mcy - 2ggx - ghy - 3cxx - 2gxy - hyy = 0, 

3{dg)y + Scdx - 2hhy - hgx - 3dyy - 2hxy - = 0. (156) 

We will use these criteria to determine the linearizability of the scalar third order semi- 
linear ode. 

3 Derivation of the linearizability criteria for the third order 
ode 

On differentiating (13) for the single differential equation, writing the independent vari- 
able as X and the dependent variable as y, we get 

y'" + (3cy'2 _ 2gy' + h)y" + Cyv'^ + (c^ - gy)y'^ - {gx - hy)y'^ + {hx - dy)y' - 4 = 0, (16) 

which is a total derivative. We could retain the equation in this form or use the original 
to remove the second order term and write the equation as quintically non-linear in the 
first derivative: 

y'" - 3c2y'5 + {beg + ey)y'^ - {4eh + 2^^ + - Cx)y'^ 
-{2gd + -hx + dy)y'^ - {2dg + + dy - hx)y' + {dh - dx) = 0. (17) 
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As we show shortly, for the linearizability criteria it makes no difference which form we use. 
However the second, quintic, form is no longer a total derivative of the original equation 
except in the linear case. For purposes of comparison with the forms in the literature one 
needs (16) but for comparison with the previous work on the cubically semi-linear systems 
of odes (17) may be more convenient. 

For (17) the general form of the equation is 

2/"' - ay'S + (5y'^ - 7^'^ + 5y'^ -ey' + <t> = Q. (18) 

In this case for compatibility of (17) and (18) the conditions are 

a = 3c2, (19) 

I3 = bcg + Cy, (20) 

-i = Ach + 2g'^ +gy-c^, (21) 

5 = Zed + Zgh + hy — g^, (22) 

e = 2dg + h^ + dy- h^, (23) 

(f) = dh-d^. (24) 

On inversion, the first four equations provide definitions for the four coefficients in (13), 
(14) and (15b) for n = 2, and the next two provide consistency constraint conditions: 

c = V^/VS, (25) 

g = {f3- cy)/5c, (c 7^ 0), (26) 

h = {-f-2g^ -gy + c^)/Ac,{c^O), (27) 

d={6-3gh-hy + g^)/3c, (c / 0), (28) 

€ = 2dg + h^ + dy-h^,{c^Q), (29) 

(t) = dh- d^, (c ^ 0). (30) 

In the case c = 0, clearly a = P = 0, and the equation becomes cubically non-linear. Now 
there can be different choices of g for a given 7, h for a given 5 and choice of g, and d 
for a given e and choices of g and h. There is then only one consistency condition for the 
various choices. An example will be given in the next section to illustrate this case. 

Hence we have the following theorem. 
Theorem 1: Equation (18) is linearizable if the linearizability criteria (15) are satisfied, 
where the coefficients are given by (25) - (30) (c ^ 0), with h = a — 2e, g = f — 2b. 

The general form corresponding to (16) is 

y'" + (^22/" - Aiy' + Ao)y" + B^y'^ - B^y'^ + - B^y' + Bo = 0. (31) 
The identification of most of the coefficients is easier in this form. We have 



c = A2/3,g = Ai/2,h = Ao, 



(32) 
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but d is not obtained from here. The constraint conditions arising are: 

B4 = A2y/3, B3 = Aiy/2 - B2 = Aoy - Ai^/2, (33) 

which are also easier to check than the corresponding equations in the other form. There 
are two differential equations for d that yield its value up to an arbitrary constant 

d = - J B2dx + k{y) = j{Bi- Aox)dy + l{x). (34) 

The constant can be determined by requiring consistency of the Christoffel symbol F^^. 
This is why it becomes more difficult to compute with this form of the equation. Thus we 
have the following theorem. 

Theorem 2: Equation (31) is linearizable if the linearizability criteria (15) are satisfied, 
where the coefficients are given hy (32) - (34) (c^ 0), with h = a — 2e,g = f — 2b, after 
requiring consistency of the Christoffel symbols with the deduced metric coefficients. 

The four linearizability conditions (15) are stated in terms of the 6 coefficients a, / 
and not the four coefficients c, d, g, h. Thus there is degeneracy in the choices available. 
Any choices of a and e for a given h, or / and b for a given g, compatible with the metric 
coefficient relations (12) are permissible. For each such choice we would get corresponding 
linearizability conditions. 

For example, assume that b = e = 0, then 

Px = -2{ap + dq),py = 0,q^ = -{aq + dr), (35) 

% = -{cp + fq),rx = 0,ry = -2{cq + fr). (36) 
This is a consistent set of requirements. The other four conditions become 

hy = cd, Cx = ch, dy = dg, gx = hy. (37) 

The first two of these are 

ihy = 6- 3gh + gx, 
ax/2a = [50V3a^7 + 25a^/2 + V3{-2V3/3 + Uyf 
+5 V3(-a^ + aayy + 2V3a{Pay - al3y))]200a^/'^ . (38) 
Hence the first and last of the four conditions are 

gx = hy = 5/3 - gh. (39) 

Since h itself depends on the second derivative of a relative to y, the last condition above 
is a third order derivative constraint. 

Alternatively, choosing a = / = 0, we would get 

Px = -2dq,py = -2{bp + eq),qx = -{hp + eq + dr), 

1y = -icp + bq + er),rx = -2{bq + er),ry = -2cq, (40) 

and 

gx = hy = 2cd - gh/2, 2c^ + gy = -2(6^ + ce), 2^ + gx = -{dg + h^/2). (41) 
These are clearly more complicated and will not be used. 
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4 Examples 

We present some examples of third order equations that can be linearized by our procedure. 

1. Here we choose e = / = so that g = —26, h = a. Further, let us try to obtain a 
lower degree of the equation in the form that there is no second derivative. As pointed out 
earlier this requires that c = 0. To construct this example we choose the metric coefficients 
p = A{x)y'^,q = B'{x)y, r = 2B{x). Writing pr - = Ay'^ and + A'B'/2 - AB" = A, 
we get h = a = —{ln^fK)',g = —2b = 2/y,c = 0,d = Ay/A,e = / = 0, where the prime 
here refers to differentiation with respect to only x and not relative to both x and y. Notice 
that the degeneracy of the case when c = has been taken care of by selecting the metric 
coefficients appropriately. Here, choosing A = c\e~^^,B = C2e~'^^ the equation becomes 

/' - 6y'Vy' + S^y'Vy - (fc' - 5/)y' + kly = 0, (42) 

where ci = 2lc2- This is amenable to reduction by the method of [14] and is not of a form 
given in [16] or [7]. 

2. Here we construct the example by differentiating the linearizable equation 

y" + xy'^ + -y' = 0, (43) 

X 

which yields 

y'" _ 3^2^/5 _ 7^/3 _ ^ (44) 

This has c = x,d = 0, 3 = / — 26 = and h = a — 2e = 2/x. We choose 6 = 
and e = — 1/x. Then / = and a = 0. It is easily verified that the conditions are 
met and that it is not in any of the classes of [16], [7] or [14] and cannot be linearized 
according to the methods therein. It can easily be reduced to the simplest linear form by 
the transformation u = x cos y,v = x sin y. Here the solution is in terms of two arbitrary 
constants. 

3. In this example the equation taken is 

y'" - Sx^y'V/ - 3a;y'Vy' + 6y'Vy' + %'V^2/ - ^v'/^^ = 0- (45) 

Since the equation is non-trivial in that it is not in the class of Neut and Petiot [16], 
Ibragimov and Meleshko [7] or Meleshko [14], is not a total derivative, nor can it be solved 
by other simple means, we give the steps in some detail. From (25) - (30), we sec that 
c = —x/y'^,d = 0,h = 2/x,g = 1/y. Choosing a = 6 = we have / = g', e = —/i/2. The 
metric coefficients are then given by (12) to be 

p = y\l + y-^)/2, q = xy{l - y-^)/2, r = x\l + y-^)/2. (46) 

Writing the Cartesian coordinates as {u,v), the coordinate transformations are given by 
([11], eqs. (32) - (34)) 

+ = y'^ + y~'^,UxUy + Va:Vy = x{y-y~^),ul + Vy = x^(y^ -|- y"^). (47) 

This set of equations is easily solved by setting 

Ux = y,Uy = X, Vx = y'^,Vy = -xy~'^, (48) 
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which yields 



u = xy, V = xy 



-1 



(49) 



which are the hnearizing transformations. Thus the solution is 



Axy + Bx/y = 1, 



(50) 



where A and B are arbitrary constant real numbers. Notice that this equation is quintic 
in the first derivative and does depend on x,y,y', which is not in the classes of [16], [7] or 



5 Concluding remarks 

We have shown that the criteria for linearizability of a third order semi-linear ode that 
were available in the literature [16, 7] exclude classes of third order equations that can 
be obtained by differentiating second order linearizablc equations. By differentiating the 
general second order cubically semi-linear ode we can check the linearizability of the class 
of third order semi-linear equations that are products of a quadratic factor in the first 
order times the second order term and quartic in the first order terms. The class can be 
more conveniently expressed as a quintic in the first order. In this case it is no longer a 
total derivative of the original second-order equation (except in the linear case) so that it 
is now hidden that this third order equation arises from a second order one. As such, our 
procedure gives a non- classical solution. The criteria are explicit for the quintic form but 
become degenerate in the case that the factor reduces to a cubic. It is worth noting that 
this proccdTirc docs not allow for a quartic factor only. Some examples were given. 

The procedure presented here leaves freedom of choice of the coefficients of the system of 
quadratic equations that yields the cubic. The reason is that there were six coefficients of 
the two-dimensional quadratic system, while only four combinations of tlicm enter into the 
third order equation. There are, indeed, six coefficients in the general third order equation 
of the desired type. The extra two coefficients yield additional constraint equations for 
the compatibility of the cubic form with the derivative of the quadratic form. We then 
only have to check that the Lie criteria are satisfied by the corresponding second order 
equation. 

Since [16] and [7] used the classic Lie procedure and obtained all the classes possible 
for the third order ode to be linearizable by point transformations, the question arises how 
we have got additional classes. The answer is that we are not using point transformations 
directly. We use point transformations and then differentiate. This could be something 
like contact transformations (as the derivative is also involved) but it is not quite a contact 
transformation procedure either. Note that our procedure does not guarantee that there 
will be three arbitrary constants. We only guarantee two. The procedure of [16] does 
guarantee three. It would be interesting if a procedure could be found that involves only 
a single constant. 

Our procedure is also interesting in that it can yield criteria for systems of third order 
odes. Work is in progress for a system of third order cubically semi-linear system [13]. 
The equivalent of the present work for a system is in principle possible but appears too 
messy at present. Better algebraic computing software may make it more manageable. 



[14]. 
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Another line of work that will be useful is to find methods for the degeneracy in our 
procedure to be removed. Alternatively, we need to determine all classes that allow for 
the degeneracy, where it cannot be removed. 
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